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Abstract 

A bound state problem in a topologically massive quantum electrodynamics is 
investigated by using a non-perturbative method. We formulate the Bethe-Salpeter 
equation for scalar bound states composed of massive fermion and anti-fermion 
pair under the lowest ladder approximation. In a large mass expansion for the 
(anti-) fermion, we derive the Schrodinger equation and solve it by a numerical 
method. The energy eigenvalues of bound states are evaluated for various values 
of a topological mass and also a fermion mass. Then we find a novel logarithmic 
scaling behaviour of the binding energy in varying the topological mass, fermion 
mass and also a quantum number. There exists a critical value of the topological 
mass, beyond which the bound states disappear. As the topological mass decreases, 
the energy eigenvalues of the bound states, which are negative, also decrease with 
a logarithmic dependence on the topological mass. A Chern-Simons term gives the 
bound system a repulsive effect. 
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1 Introduction 



A peculiar type of quantum field theory exists in (2+l)-dimensional (generally odd- 
dimensional) space-time. It is the theory that has a Chern-Simons term in the action [[[]. . 
While the Chern-Simons term may be added to the action by hand, the term can be in- 
duced in an effective action by a vacuum polarization effect of a fermion as the anomaly [[J . 
Our understanding on the theory has progressed by studying the topological structure due 
to the Chern-Simons term [f| 

In addition to the theoretical interests, two important phenomena were discovered in 
planar condensed matter systems. One is the quantum Hall effect [f| and the other is 
the high-T c superconductivity ||. Both phenomena may be thought of as realizations of 
the macroscopic quantum effects in the planar systems in which electrons have a strong 
correlation. It is plausible that these phenomena may have their origins in the dimensions 
of the space-times. It is challenging to ask whether and how the Chern-Simons term is 
related to these phenomena discovered in the two-dimensional electron systems. 

These situations have activated the study of the (2+l)-dimensional quantum field the- 
ories. In fact, there appeared many approaches to understand these macroscopic quantum 
effects by using the (2+f )-dimensional quantum field theories and by paying special at- 
tention to the effects of the Chern-Simons term. For example, the (2+f )-dimensional 
quantum electrodynamics (QED 3 ) was used to explain the quantum Hall effect [0] and 
QED 3 with the Chern-Simons term (Maxwell-Chern-Simons QED 3 ) provided an anyon 
model H which was expected to give us an essential mechanism for the high-Tc super- 
conductivity. These investigations produced important results and are still in progress. 

What is the physical meaning of the Chern-Simons term? In the (2+f )-dimensional 
gauge theories, the behaviors in the infrared regions seem to be unstable at least in the 
perturbative treatment. The Chern-Simons term gives the gauge field a topological mass 
without breaking the gauge symmetry so that the term may rescue the theory from the 
infrared catastrophe. This is an original motivation of including the Chern-Simons term 
in the action. || 

Reminding those studies, we think that one of the most important problems which are 
still unclear is what is a role of the Chern-Simons term in a nonperturbative region. It is 
our purpose to clarify how the Chern-Simons term affects the nonperturbative dynamics. 
In previous works |ITJ [IT], |T2|, |TB|, II |, a dynamical mass generation of a fermion has been 



investigated by using the Schwinger-Dyson equation. It is an important phenomenon 
induced by a nonperturbative effect in one-body system. In this paper, we study a bound 
state problem which is another important nonperturbative phenomenon appeared in two- 
body system. We are interested in how the Chern-Simons term affects the bound system 
through the nonperturbative dynamics. 

In Ref. [f[I|| , a fermion-anti-fermion bound system in QED 3 was studied. In the model 
used there, the (anti-) fermion had four components and the gauge field did not have the 
topological mass. On the other hand, a binding problem of fermion- fermion system |f6| 



f 



was investigated in Ref. [I7|. A two- component fermion was used and the Chern-Simons 



term was included in the action of gauge field. The effective fermion-fermion low-energy 
potential was evaluated in the lowest order perturbation. However, the result was claimed 



to be incorrect [IS, IS 



As mentioned above, what we are interested in is the effect of the Chern-Simons term 
in the nonperturbative dynamics. It is still unclear and controversial. We think that it 
is important to perform the microscopic nonperturbative calculation based on the first 
principle. For this purpose, the best theoretical tool is the Bethe-Salpeter equation which 
makes it possible to perform the systematic nonperturbative analysis of the structure 
of the bound system. As a system suitable for our purpose, we use the two-component 
massive fermion and anti-fermion in QED 3 with the Chern-Simons term. This model has 
the topological nature and also has applications in the condensed matter physics. 

This paper is organized as follows: In Sec. 2, we explain the quantum electrodynam- 
ics with the Chern-Simons term. The Bethe-Salpeter equation is formulated under the 
lowest ladder approximation in Sec. 3. In Sec. 4, we derive the Schrodinger equation 
using the large fermion mass expansion. After that, by solving the Schrodinger equation 
numerically, we evaluate the binding energy for the bound states and we find a specific 
logarithmic scaling behaviour of the binding energy in Sec. 5. Finally, we conclude on 
the results with discussions in Sec. 6. 



2 Maxwell- Chern-Simons QED3 

The Lagrangian density of QED3 with the Chern-Simons term, which is called the 
Maxwell-Chern-Simons QED3 or the topologically massive QED3, is expressed as 

L = -~F^F^ + ^ p A^A p - ^{d^f + 4,{i ? - e / - m)iP , (1) 
4 z la 

where e is the coupling constant and a is the gauge-fixing parameter, ip is a two com- 
ponent fermion field with an explicit mass m, which belongs to the irreducible spinor 
representation in (2+l)-dimensions. In the dimensions, the square of the coupling con- 
stant e 2 has a dimension of mass. The second term on the right-hand side of Eq. ([!]) is 
the so-called Chern-Simons term. This term has a topological meaning as a secondary 
characteristic class J|, |(| and provides the topological mass \i to the gauge field without 
breaking the gauge symmetry. The free propagator of the gauge field iD pv {p) derived 
from Eq. (|l|) is 

■n f \ ■ 1 ( PnP»\ . 11 P ■ VvPv 

iD^(p) = -1— — + fi— -— e^ p p p - ia^- . (2) 

p z — fi z \ P J P ~ P> p p 

We find a massive pole at p 2 = /1 2 in the propagator so that /i can be regarded as a mass 
of the gauge field. 
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The Dirac matrices in (2+l)-dimensions are defined by 

7° = o"3, 7 1 = MTi, 7 2 = ^2 , (3) 

with diag((7^) = (1, —1, —1). oVs (i = 1,2,3) are the Pauli matrices. In this representa- 
tion, there is no matrix which anti-commutes with all of the 7-matrices so that we cannot 
define the chiral transformation. Instead, the parity symmetry [111 3 forbids the mass 
terms of both the fermion and the gauge field. The Lagrangian density ([!]) breaks the 
parity symmetry in both the gauge and fermion sectors. 



3 Bethe-Salpeter Equation 

The Bethe-Salpeter amplitude XB(x a ,x b ) is defined as 

X %{x ai x b ) = ^\TW{x a W(x b ))\k) , (4) 

where \k) is the fermion-anti- fermion bound state with the total momentum and mass 
M. The mass shell condition is k 2 = k\ — k 2 = M 2 with < M < 2m. The integral 
Bethe-Salpeter equation |2l], p2| for XB(x a ,x b ) is given by 

XB(x a ,x b ) = -e 2 J d s z a d 3 z b iS' F (x a - z a )YxB(z a , z b )Y 

xiS' F (z b - x b )iD'^{z a - z b ) , (5) 

where iS' F and iD'^ v are full propagators of the fermion and gauge fields respectively. 

Now we limit ourselves to solve Eq. (||) in the lowest ladder approximation. Then iS' F 
is replaced by the free fermion propagator iS F as 

iS F (x a -x b ) = if 4^ e ~ tk(Xa ~ Xb) Mk), ^S F (k) = i J ft+ ^_ . , (6) 
and iD'^ by the free gauge field propagator iD^ as 

iD^e-^-^, (7) 

where iD^ u (k) is given by Eq. (Q). 

The dependence of XB(x a , x b ) on the coordinates of the centre of mass can be factored 
out as 

X l i(x a ,x b )=e- lPX X %(x,P), (8) 
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where X = (x a +x i )/2 is the coordinates of the centre of mass and x = x a — Xb the relative 
coordinates. P is the momentum of the centre of mass. By the Fourier transformation of 
the relative coordinates x, we obtain 



XB^ X ai x b) 



-IPX 



d s p 



(9) 



where p is the relative momentum. 

We may rewrite the Bethe-Salpeter equation by using an amputated Bethe-Salpeter 
amplitude T which is defined as 



iSf 



r 



+ 



-i 



Xb(p,P) 



(10) 



From Eqs. @ ~ ([T0|), we obtain 
d 3 p 



(2tt; 



3 iD^(q-p)^iS F (^-+ j)T(p,P)iS F (-£+ th v 



'111 



Now the problem is to find the amputated Bethe-Salpeter amplitude V by solving Eq. (|TT|). 
Because of the Lorentz invariance, the amplitude which is a bi-spinor, can be decomposed 

as 



T(q, P) = M{q, P) + N(q, P) 4 + F(q, P) f + G(q, P)e^%P ulp , 



(12) 



where M, N, F and G are the functions of the Lorentz scalar variables q = y/q^q 11 , P = 

J PfiP^ and Pq = P^. The equations for the four unknown functions M,N,F and G 
can be obtained by taking the trace of the function V as 



M(q,P) 
N(q,P) 

F(q,P) 

G(q,P) 



^r[T(q,P)] , 

1 PV - (P • q)P a 

2 P 2 q 2 - (P ■ q) 2 

1 q 2 P a - (P • q)q a 

2 P 2 q 2 - (P • q) 2 
1 e^g^P/3 



tr[ 7(7 r(g,P)] 
tr[ T(j r(g,P)] 



2 P 2 g 2 - (P • g) : 



tr[ 7ff r( 9 ,P)] . 



(13) 
(14) 

(15) 

(16) 



In the following section, we rewrite the obtained four coupled integral equations into 
the non-relativistic form using the large fermion mass expansion. This leads to the 
Schrodinger equation which will be used actually to study the bound states. The non- 
relativistic treatment reduces the complexity of the equations as we will see in the next 
section. 
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4 Schrodinger Equation 



We consider the bound state in the rest frame where the momentum of the centre of mass 
is 

P M = (2m + E B ,0,0). (17) 

E B indicates the binding energy and m the (anti-) fermion mass. If the system makes 
bound states, E B should be negative because of a stability. We replace all P M 's in Eq. 
(|TT|) by Eq. (|17|). Then the fermion propagator in the right-hand side of Eq. (|TT| ) can be 
approximated as 

±l+p' + m ^ 1±7 o 



(±f +p) 2 -m 2 + ie E B ±2 Po - ±f + ie' 

in the large limit of the fermion mass, pP/m 2 term must be kept in the equation since 
this term may play an important role for the convergence of the integration over p in the 
right-hand side of Eq. (ITT1). The amputated Bethe-Salpeter amplitude T in the rest frame 
of bound states can be written as 

T(p,P) = M(p, P) + N(p, P) p 1 + F(p, P)(2m + E B )^° 

+G(p,P)e^p l (2m + E B )^ , (19) 

by substituting Eq. ( |T7| ) into Eq. (|T2"D. Using Eqs. (|T8| ) and (|T9|), we find the non- 
relativistic Bethe-Salpeter equation for the amputated Bethe-Salpeter amplitude T as 

T(q,P) = 2e 2 J -^iD^q-p^il+l^YK^P), (20) 

where 

KfaP) = {N(p,P)p J + (2m + E B )G(p,P)e i03 p i } 1 ^ 1 ^, (21) 

K(p,P) = (E B -Lf + ief-Apl. (22) 

Now, we can derive the four coupled integral equations for four unknown functions M, N, F 
and G from Eq. ©. 

After tedious calculations^ we find that the essential coupled equations are those 
of iV and G. M and F are calculated from iV and G. Furthermore, we find that the 
independent unknown function is only T(q, P) defined by 

T(q, P) = N(q, P) - i(2m + E B )G(q, P) , (23) 



1 The details will be appeared elsewhere. [g3| 
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which is required to satisfy 

T(9 ' p) = / Wf R{q - p) ^-P-^^^§^-y (24) 

where 

The symbol "*" was defined as q * p = e^q^Pj. Thus by taking the non-relativistic limit, 
the number of the independent equations reduces from four to one. Next we derive the 
Schrodinger equation from the non-relativistic Bethe-Salpeter equation (|24|). Using the 
residue theorem, we perform the integration over p° in Eq. fl24l) . We consider poles in the 
function K(p, P) defined by Eq. fl22|) and get 

T(q ,q,P)= 1 - [ ^- 2 R(q - p* Q~ p)(~g " P + iq* P) ^ f P) , (26) 

where = (E B — p 2 /m)/2. If we fix g to = (i? 5 — q 2 /m)/2 in the both sides of 
Eq. (|26|) and also neglect q 2 /2m and p 2 /2m for the large fermion mass m again, R(q — p) 
reduces to a simple form as 

Then, as the self-consistent equation, we obtain 

T{q 0l q,P) = ~e jtt^^J^ — — 2 # r- 28 

7 (27r) 2 g 2 (q - p) 2 + fi 2 E B - ^ + «e 

Now we have to define the Schodinger wave function ip B (p) of the bound states from 
T(p B ,p,T). It is natural to assume that there exists an effective theory for the bound 
states. The Hamiltonian of the effective theory, from which the Schrodinger equation is 
derived by a variational principle, should be hermite since an energy eigenvalue has to be 
real. The hermiticity decides that the wave function should be defined as 

* {P) ~ (E B -^ + ie) ' { } 

up to a multiplicative constant .0 Finally we have 

(E B ~ ~ + ie)V^(?) = -e 2 I 7^t4 ^ ^f. ^ *(p) ■ (30) 



m ' J (2ir) 2 \q\\p\ (q - p) 2 + fi 2 

Notice that the term with ^-product drops in Eq. (|30"D because the wave function depends 
on only \p\ . This is the Schrodinger equation in the momentum space which we solve in 
Sec. 5 to obtain the binding energy . 

2 See Ref. S for the details 
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5 Numerical Analysis 



5.1 Control parameter 

In the Maxwell-Chern-Simons QED3, the square of the coupling constant e 2 has the 
dimension of mass. The theory is under a control of dimensionless parameters defined by 

V = ~v m "- ( 31 ) 

These two parameters, not three, are independent in the theory Further more, by defining 
the dimensionless momenta q and p, and dimensionless binding energy Eb as 



iB 



we can rewrite the Schrodinger equation (|30|) to 



/ f?B ^ 2 V? . ■ wB,^ 1 f d2 P 1 Q-P ,B,*s /qq\ 



We can see that only the combination of parameters as 



^ ^ (34) 



me 2 m 



fully controls Eq. (|33|). Namely, if one chooses a certain value of /1 2 / (me 2 ), the dimension- 
less energy _E S and the wave function ijj B (q) are decided. To obtain the eigenenergy E B 
from the dimensionless E B , we need to fix one more parameter, m or fx. In the followings, 
we draw the results in pictures by fixing m to see the topological mass dependence of E B 
and also by fixing \x to see the fermion mass dependence. 



5.2 Topological mass dependence 

What we are most interested in is to find how the topological mass affects a property of 
the bound system. It can be seen in the shift of the energy eigenvalue by varying the 
topological mass. Especially, it is very important to clarify whether the effect by the 
topological mass is attractive or repulsive and how strong the effect is. 

In Fig. 1, we show the topological mass dependence of the dimensionless binding 
energy E B /e 2 . The vertical axis indicates the binding energy and the horizontal axis is 
the topological mass fi, in the logarithmic scale. As the topological mass jl increases, the 
absolute value of the binding energy becomes smaller and the number of the bound states 
decreases. We find that there exists a critical value of the topological mass fx where the 
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Fig. 1. The dependence of the dimensionless binding eigenenergy E B je 2 on the dimen- 
sionless topological mass fi with the dimensionless fermion mass rh = 5.0 x 10 4 . Open 
rhombuses indicate the calculated binding energies E B /e 2 . The dotted lines are straight 
with the slope 1/(2-k) as given by Eq. (pop. 



bound states disappear, for each a quantum number. For example, the critical value of 
topological mass fi is numerically 44.48 — 44.49 in the case of the ground state. 

On the other hand, we can see another interesting feature of the topological mass 
dependence of the binding energy from this result. In the region where the binding energy 
is lower than about —0.3, the topological mass dependence of the binding energy is almost 
linear in the semi-logarithmic figure and the value of its slope is l/(27r) at least for lower 
energy states. Namely, the dependence of the binding energy seems to be approximated 
by a function form as 



e 2 2tt 



■2 o > ( 35 ) 



where the index n is a quantum number indicating the energy level. is determined 
so as to reproduce the calculated value of E B /e 2 at p, = 0.4 for each n in Fig. 1. Before 
we discuss about the meaning of this dependence in Eq. ([35]), we show the fermion mass 
dependence of the binding energy. After that, we shall return to this point. 

5.3 Fermion mass dependence 

In Fig. 2, we show the fermion mass dependence of the binding energy. The vertical 
axis indicates the dimensionless binding energy E B /e 2 and the horizontal axis shows the 
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fermion mass rh in the logarithmic scale. The topological mass p, is fixed to be 1.0. As we 
can see from this figure, the absolute value of the binding energy increases as the fermion 
mass m increases. We find the interesting feature again from this figure that the results 
are well approximated by the function as we have seen in the case of the topological mass 
dependence. In the region of the binding energy lower than about —0.3, the fermion 
mass dependence is almost linear with the slope l/(47r) in the semi-logarithmic figure. 
Therefore, the fermion mass dependence of the binding energy seems to be approximated 
by the function form as 



e 2 4tt VeV n ' 



(36) 



where is determined by the value of E B /e 2 at m/e 2 = 5.0 x 10 5 for each n in Fig. 2. 



fermion mass dependence (topological mass=1 .0) 
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Fig. 2. The dependence of the dimensionless binding energy E B /e 2 on the dimensionless 
fermion mass rh with the dimensionless topological mass fi = 1.0. Open rhombuses 
indicate the calculated binding energies E B /e 2 . The results can be fitted by the dotted 
lines defined in Eq. (|3T 



5.4 Logarithmic scaling of binding energy 

Here, we return to the topological mass dependence of the binding energy expressed by 
Eq. (|35D . By combining Eqs. (|35| ) and (|36|) , we find that the dependence of the binding 
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energy on m and ft can be unified to the expression as 



E 

e 



Air 



In 



me' 



(37) 



where is a constant to be consistent with Eqs. (|35|) and (|36|) . We can see that 
the dimensionless binding energy E B /e 2 is completely determined by the value of the 
parameter fi 2 / (me 2 ) which is only one essential parameter included in the theory as seen 
in Sec. 5.1. 



quantum number dependence (fermion mass m=5.0e4) 




quantum number n 

Fig. 3. The dependence of the dimensionless binding energy E B /e 2 on the quantum 
number n. The results can be fitted by the dotted lines defined in Eq. (|38D with a = 
2.1157 ~ 2.1505. 



In addition, we see the dependence of the binding energy on the quantum number n in 
Fig. 3. We can fit the data points for lower energies by a logarithmic function. Combining 
the fitting with Eq. ([37]), we obtain 




where a is a fitting parameter used in Fig. 3. Compared with the fittings in Figs. 1 and 
2, the fitting in Fig. 3 results us a deviation of the values of a for each fitting line. For 
n — 1 ~ 5, the fitting parameter varies as a = 2.1157 ~ 2.1505. 

The most interesting point of Eq. (|3"8|) is that E B /e 2 is scaled by the logarithmic 



function of n /{me ) and also n. If we make a shift as /(me ) — > Xfi /(me 



or n 
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An where A is a constant, then E^/e 2 is changed as E^/e 2 — > E^/e 2 + l/(47r)lnA or 
En/ e2 ~^ En/ e2 + l/(o7r)lnA. Thus the logarithmic scaling behaviour means that the 
multiplicative shift of the parameter or the quantum number results the additive shift of 
the binding energy. 

The scaling has been found only through the numerical analysis in this paper. We 
may expect any simple explanation about the approximated logarithmic dependence of 
E B on the unique parameter /i 2 / (me 2 ) and also the energy level n. The coefficient l/(47r) 
is also suggestive for us since it might be a signal of a topological nature in the theory. 
These will be studied in future investigations. 

6 Conclusion 

We have investigated the bound states of the fermion-anti-fermion system in the three- 
dimensional QED with the Chern-Simons term. We used the two-component fermion 
field which belongs to the irreducible spinor representation in (2+l)-dimensions. We 
have formulated the Bethe-Salpeter equation which makes it possible to perform the 
systematic non-perturbative evaluation based on the first principle. We have performed 
the calculations in the non-relativistic limit constructing the Schrodinger equation in the 
momentum space under the large fermion mass expansion. By solving this Schrodinger 
equation numerically we investigated how the Chern-Simons term affects to the fermion- 
anti-fermion bound states. 

We have discovered some important features of the effects of the Chern-Simons term in 
the numerical analysis. First, we have found that the absolute value of the binding energy 
becomes smaller as the topological mass /x larger. This fact indicates that the effect of the 
Chern-Simons term is to weaken the binding of the fermion-anti-fermion bound states. 
We have found that there exists a critical value of the topological mass where the bound 
state disappears. 

Secondly, we have found the interesting features of the binding energy that in the 
region where the energy eigenvalue is smaller than about —0.3, the dependence of the 
binding energy on the topological mass /i, the fermion mass m and the energy level 
n is almost linear in the semi- logarithmic plots. We find that this dependence can be 
expressed as E^ = e 2 /(An) log(/i 2 /(me 2 ) x n 4 / a ) with a = 2.1157 ~ 2.1505 for the energy 
levels n — 1 ~ 5. Thus the binding energy has a logarithmic scaling behaviour for varying 
the topological mass, fermion mass and also the quantum number. 

In future investigations, first we should make clear the meanings of the logarithmic 
scaling, the coefficient l/(47r) and the critical value of the topological mass. We should try 
to obtain any analytical information which helps us to understand the meanings of these 
behaviors of the numerical results. Further, we can investigate this problem using the 
relativistic formulation itself without making the non-relativistic approximation. We will 
develop the combined treatment of the Schwinger- Dyson equation and the Bethe-Salpeter 
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equation in order to incorporate the effect of the dynamical mass generation. In addition, 
we think that one of the most interesting subjects which should be done in future is to 
extend this study to the bound states of the fermion-fermion system |16], |24|. This is 
important due to the possible connection to the high-Tc superconductivity observed in 
the condensed matter systems. 
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